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Abstract 

A hybrid model which allows to interpolate between the (original) Regge approach and 
dynamical triangulation is introduced. The gained flexibility in the measure is exploited to 
study dynamical triangulation in a fixed geometry. Our numerical results support KPZ 
exponents. A critical assessment concerning the apparent lack of gravitational effects in 

two dimensions follows. 
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1 Introduction 



Two numerical methods, both derived from Regge calculus, have been used to investi- 
gate non-perturbative quantum gravity: dynamical triangulation (DT) for fixed, equal link 
lengths, called standard dynamical triangulation (SDT) henceforth, and the (original) Regge 
approach (RA), varying the link lengths of a simplicial lattice with fixed incidence matrix. 
The latter suffers from the same fundamental problem as continuum formulations: the cor- 
rect measure is unknown. On the other hand, the SDT approach is claimed to be free of 
those troubles since the summation of triangulations is thought to provide the unique, correct 
measure of geometries. In the following we focus on 2D Ising spins coupled to the geometry, 
as in this case a number of exact results are known. 

For 2D Ising spins the SDT method is known to be equivalent to exactly solvable matrix 
models [1, 2]. The fact that KPZ theory [3] reproduces precisely the the same critical 
exponents is widely regarded as a "proof [4] that the SDT method is correct (not only in 
2D). In contrast the RA has failed to obtain the KPZ results [5, 6]. Instead consistency with 
Onsager's exponents for the flat space Ising model was found. Recently, the scenario has 
become further diffused by the claim of Vekic et al. that KPZ exponents are observed for a 
flat space Ising spin system, where variable incidence matrices are used [7]. 

Let us introduce a hybrid model which allows to interpolate between RA and SDT. The 
path integral reads 

Z No (P,\,a) = J2W T j D{l T }J2e~ E (la) 

with 

£ = 2( Mf + a (y) + ^_ 4 , T (fa-i_f«) 2 (14) 

Here s 8 = ±1, i = 1, N are the usual Ising spins, and i,j label sites of the lattice. The 
first sum is over distinct triangulations T of a given topology, weighted with Wj > 0. The 
adjacency matrix is defined by 

( 1, if i and j are neighbours, 
13 I 0, otherwise. 
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Now l T labels the links / 8J for which the adjacency matrix T^j^ is non-zero, where and 
are the sites connected by /. We use t T to label the triangles of the triangulation T. The 
area Aj , associated with a site i is taken to be baricentric 

where A t r denotes the area of the triangle t T . By n 8 we denote the number of triangles 
encountered in the sum, i.e. the coordination number. The deficit angle at site i is defined 
as 

Sf = 2vr - £ 9i{f), 

where 9i(t T ) is the dihedral angle associated with the triangle t T . Similarly, the area A/t, 
associated with the link l T , is defined as 

a i t = \ E A t T - 

6 t T Jl T 

Finally, we use Ni and N 2 to denote the numbers of links and triangles, respectively. Choos- 
ing Wt and D{l T } to interpolate between RA and SDT, we expect that the model can shed 
new light on the reasons for the occurance of different universality classes, characterized by 
Onsager versus KPZ exponents. For the RA all weight factors but one are zero, whereas for 
DT all weight factors are identical, for instance Wt = 1. In this paper we focus on DT with 
the measure 

D{F} = Y[T ij dl ij 6{U j -c ij ), (2) 

ij 

where (c 8J ) is a N X N matrix of constant link lengths. With 

Cij = 1 (3a) 

the SDT approach is obtained. Another notable choice is 

Cf, = geodesic distances between vertices i,j in a fixed geometry. (36) 

With (3b) it is easy to see that the geometry (of the starting configuration) stays fixed. Vari- 
ations due to changes in the incidence matrix Tf, are compensated by appropriate changes of 



the link lengths. Clearly, the general measure in (la) allows to interpolate smoothly between 
(3a) and (3b), for instance by giving the two length assignments different weights which add 
to one. Such a parametrization allows to study the onset of gravitational dressing when 
proceeding from (3b) to (3a). Correspondingly, one may expect the critical exponents to 
vary between Onsager for (3b) and KPZ for (3a). However, in the next section we present 
simulation results with the measure (3b), already supporting KPZ exponents. The DT sce- 
nario becomes thus very similar to the RA results reported in [5, 6]: Turning "gravity" on 
or off seems to have no effect on the critical exponents. A critical assessment of these results 
is tried in section 3. Summary and conclusions are given in section 4. 

2 Simulations 

We present simulation results with the measure (3b). The action is given by (lb) with a = 0. 
Since we consider ensembles of fixed area the parameter A is of no significance. As in previous 
simulations [8, 9, 10, 11] with the 2D DT approach we rely on the spherical topology. 

2.1 Algorithm 

We first describe how both Tf, and c 8J can be generated through a Metropolis algorithm. We 
start with a geometry specified by a triangulation T° with equilateral triangles: 

T? 3 = 1 c t3 = 1. (4) 

All other triangulations are generated as for SDT by applying successively link-flip moves. 
In difference to SDT we change the lengths of the flipped links to get c 8J as defined by (3b). 
The procedure is as follows: The two triangles affected by the link flip are embedded in R 2 
as shown in figure 1. The vertices are labelled as z, j, k } I. The length c 8J of the flipped link 
is calculated from the known lengths, c 8J = Cij(cik,cn,Cjk,Cji), using the Euclidean distance 
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in flat space 



c li c ii + c ik , „ 

= ' h = ^ * k ~ 1 

-2 = — , v, = ^c) k -A 



'(yi + 1J2) 2 + (xi - x 2 ) 2 . 

It is easily seen, that all (geodesic) distances, the deficit angles and the (total) area are 
preserved by this flip. The geometry stays unchanged. Since link flips may successively lead 
to large lengths / 8J it is necessary to use a cutoff Iq so that only configurations with / 8J < Iq 
are allowed. Also the algorithm has to reject certain flips to enforce that the area stays 
invariant, see figure 1. In addition we reject flips which lead to triangles with area zero. 

In the next subsection we present numerical results for two different geometries. To 
construct T°, we begin with the simplest triangulation of the 2-sphere: the surface of a 
tetrahedron. This provides for a list of N = 4 vertices, Ni = 6 links and N 2 = 4 triangles 
together with incidence tables describing how they are connected. Then one divides succes- 
sively the links applying Alexander [12] moves, which increase N by one, Ni by three and 
N 2 by two. After each move the new simplices are added at the end of the incidence tables 
and the procedure is iterated until the desired number of vertices is reached. Increasing the 
index of the link which is split sequentially gives geometry Gl, while increasing it in steps 
of two results in Gl. Cutoffs lc = ^c = 5-4 are used for both geometries. For N = 250 the 
distributions of coordination numbers n are given in table 1. For larger N the distribution 
of the coordination numbers converges fast towards a fixed shape, which allows to speak of 
fixed geometries. 



2.2 Finite size scaling preliminaries 

For both geometries we have performed Metropolis simulations to estimate critical exponents 
through finite size scaling (FSS) analysis. Similar as in [6] we define pseudocritical points 
/3*(iVo), i = 1,2,3,4 as maxima of suitable physical quantities. They are: 
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1. The magnetic susceptibility 

X(P) = N (< m 2 > - < |m| > 2 ) , (5) 

where m is the lattice average m = Nq 1 J2i s i- 

2. The derivative of the absolute value of the magnetization: 

D m (f3) = ^— = < E >< \m\ > - < E\m\ > . (6) 

dp 

3. The derivative of the Binder [13] cumulant B = 1 — | < m 4 >< m 2 >~ 2 : 

□ . . dB , . ( < m 2 E > < m 4 E >\ 
D Z 3 = 1^ = (l-B)l<E> -2 + . 7 

dft \ < m 2 > < m 4 > J 

4. In addition the specific heat 

C v = /3 2 jV! (< e 2 > - < e > 2 ) with e = E/N 1 (8) 



is considered. The maximum values of these quantities are denoted by Xc(^o)j D™(N ), 
D^(No) and Cv }C (N ). Their FSS behavior determines the critical exponents 7, f3, v and a 
through the following equations: 

Xc = an + ai 2 i\# , = a 21 + a 22 N^ , (9a) 

1 a 

Df = a 31 + 032^0" , and °v,c = a 41 + a 42 iV(f . (96) 

Here we have neglected higher order corrections in N . The expected exponents, KPZ versus 
Onsager, are collected in table 2. As a < 0, the specific heat does not really define a 
pseudocritical point. In the N — > 00 limit, the pseudocritical points converge towards the 
critical point: 

P i c (N ) = p c + b i No&. (10) 

We employ the multicanonical method [14] to cover for each N a /3-range large enough to 
accomodate all pseudocritical points in one run. An additional advantage is that metasta- 
bilities may be avoided through excursions into the disordered phase. 



2.3 Numerical results 

Tables 3 and 4 give an overview of our statistics and simulation results for geometry Gl and 
Gl. The statistics is given in sweeps, where one sweep is defined by updating in the average 
1/2 of the links and 3/2 of the spins. Measurements were performed every ten sweeps, except 
for iVo = 4500 for which they were performed every twenty sweeps. Error bars are calculated 
with respect to twenty jackknife bins. They are given in parenthesis and apply to the last 
two given digits. 

Figure 2 shows our multicanonical energy distribution for the N = 4500 lattice. This 
should be compared with figure 3, where we show the functions x(fi)-> 15D m (/3) and 60CV(/3) 
on a semi-log scale (the factors are chosen such that all quantities can be depicted together). 
It is obvious from figure 3 that Xc should allow a more accurate FSS analysis than D™ . The 
maximum of the specific heat shifts with No, and stays essentially constant. Figure 4 depicts 
D B (f3). Due to the higher moments in its definition (7) substantially larger relative errors 
than in figure 3 are encountered. 

We now turn to the FSS analysis. For geometry Gl we have the larger statistics due to 
simulating an N = 4500 lattice. Nevertheless, we face the problem that three parameter fits 
to the equations (9), are unstable. Substantially larger lattices would be needed, but their 
simulation is impossible within our present CPU time limitations. Therefore, our analysis has 
to rely on two parameter fits. We assume exact KPZ or Onsager exponents, see table 2, and 
ask the question whether consistent fits are then possible. To decide it we have calculated 
the goodness-of-fit [15] Q. Under the assumption that the discrepancy between the data 
and the fit function is entirely due to statistical fluctuations of the data, Q is an uniformly 
distributed random variable in the range (0:1]. Our fits rely on lattices with N > 500, as 
the iVo = 250 lattice turned out to be too small to exhibit asymptotic (large N ) behavior. 
This means, we use six data points for geometry Gl and five data points for geometry Gl. 
The results for Q are collected in table 5. It is clear that Onsager exponents are ruled out, 
the likelihood that the encountered discrepancy is due to chance is apparently < 10~ 25 for 
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geometry Gl and < 2 • 10~ 4 for G2. Only the geometry G2 fit for l/(2z/) is consistent with 
the Onsager exponent, presumably due to the large statistical errors of the corresponding 
data. On the other hand, the data are well consistent with KPZ exponents. There is one 
Q = 0.04 value, but with six independent observables the likelihood to encounter this value 
would be about 20% , i.e. their is no statistically significant discrepancy. Figure 5 compares 
KPZ versus Onsager fits for [Xc(^Vo) — Xi(-^o)]) where Xi(-^o) is the geometry Gl KPZ fit 
for Xc(iVo). 

Assuming now the KPZ value for v, equation (10) is used to estimate the critical 0. 
Smaller lattices are successively eliminated from the fits, until an acceptable Q value is 
obtained. The thus obtained C estimates are also included in table 5. For geometry Gl 
the estimates indicate that systematic errors exceed the statistical error significantly. It is 
not entirely clear whether it would need larger lattices to reconcile these values. Possibly 
the observables, in particular the derivative of the Binder cumulant, suffer from systematic 
errors which one could improve by just enlarging the statistics. For G2 we do not find such 
an inconsistency, but there is also a problem with the Df data: the goodness-of-fit Q gets 
only acceptable when the fit is constrained to the last three lattices. Six data points were 
used for all other C fits. Taking weighted averages we get 

C = 1.123 ± 0.002 for Gl and & = 1.120 ± 0.005 for G2. (11) 

Here the error bar is taken from the Xc estimates, as other estimates are from the same 
configurations, it would be inappropriate to use them for error bar reduction. 

Data for the specific heat follow in table 6. In addition to the maximum values Cv, c 
we also give the values at the infinite volume critical points: Cv(0 c ) with C estimated by 
equation (11). The Onsager value a = would imply a logarithmic behavior, and we use 
the fit 

C v ,c = a 41 + a 42 ln(iVo) 

to compare it with the fit (9b) for KPZ exponents. The Q values for these fits are collected in 
the last two rows of table 6. From the maximum values Cv, c we cannot distinguish between 
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KPZ and Onsager. As one may have expected, both fits are consistent, presumably due 
to a dominant regular part. Somewhat surprising is the bad quality of the same fits when 
Cv{ftc) instead of Cv, c is used. Onsager becomes then inconsistent. Although some orders of 
magnitude better, the quality of the KPZ fits is not too convincing either. Most likely the 
reason is that the asymptotic CV(/3 C ) behavior will set it in only for larger lattices. Assuming 
now a = — 1, it should be noted that there is now reason for the location of the maximum to 
approach f3 c . Indeed, changing the multiplicative factor f3 2 in equation (8) to another power 
would shift the position of the maximum. This is the reason, why we also listed CV(/3 C ) and 
do not attempt to estimate f3 c by means of FSS analysis of the f3^(N ). 

Another set of two parameter fits is obtained by setting the constants an in equations 
(9) equal to zero. This gives direct estimates of the critical exponents involved, which for 
sufficiently accurate data from large enough lattices will converge to the exponents of the 
infinite volume theory. We include the N = 250 lattice, if the goodness-of-fit is already 
satisfactory, otherwise, we apply the fits to our N > 500 lattices. The results are collected 
in table 7. They are a bit off from KPZ exponents, although closer than typical results 
from SDT simulations [8, 9, 10, 11]. The fact that we can perform different two parameter 
fits, which are both consistent, explains why three parameter fits are still unstable. In the 
investigated range of lattices the constants an are still important, in particular for geometry 
G2. 

We have also started to investigate geometries on the torus. Apparantly this is hampered 
by the fact that regular initial configurations give rise to a substantial fraction of flips which 
result in degenerate triangles of area zero. In the present paper we have excluded these flips, 
but this is harmless as for the considered geometries on the sphere they are of measure zero. 
From the torus we have preliminary indications that there may be geometries which belong 
to yet different universality classes. The basic problem remains, to disentangle gravitational 
dressing (if any) from fixed geometry DT interactions. Future investigations may provide 
more insight. 
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3 Interpretation 



It is well known that the critical properties of the 2D Ising model correspond to a free fermion 
theory. It is natural to assume that (fixed geometry) DT moves will give rise to effective four 
fermion interactions. The anomalous dimension d^, where ip is is the massless interacting 
Thirring field, allows then to parametrize the critical exponents [16]: 

" = (2-<^) -1 and r l = d TH l<l. (12) 
The other exponents follow through the scaling laws 

1 - a/2 = v, a + 2/3 + 7 = 2 and 7 = (2 - -q)v, 

which in this model are easily seen to hold. The choice d^ = 4/3 does indeed reproduce 
the KPZ exponents. Our conjecture is that for a fixed geometry DT moves induce the 
appropriate four fermion interaction. Introducing then "gravitational" effects through SDT 
fluctuations of the geometry (i.e. setting the / 8J = 1 and thus A/t/(/ t ) 2 = y/Z/A) has no 
further effects on the critical exponents. 

SDT is considered a promising approach towards the quantization of gravity because 
it is conjectured [17, 18] to sample representative geometries. A major factor of support 
towards this conjecture has been the agreement with the exact KPZ exponents. As our 
results suggest that these exponents are not sensitive to fluctuations of the 2D geometry, 
the agreement between SDT and KPZ becomes trivial. Any sampling over geometries will 
do, as long as DT updates are used. A fresh look at the issue whether SDT will sample 
representative geometries becomes legitimate. 

We like to raise the concern that there are geometries which do not allow for smooth 
approximation by triangulation with equilateral triangles, even in the limit N — > 00. We 
consider as observable the deficit angle 6[p] associated with a closed path p. As is well known, 
any vector parallel transported around a closed path will in general undergo a rotation, such 
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that the initial vector and the final vector differ by a deficit angle 6[p\. (In 2D S[p] is 
the same for all kinds of vectors transported along the path p and describes directly the 
curvature of the enclosed region.) A closed path on a simplicial lattice encircles in general 
N[p] vertices and 6[p] is the sum of deficit angles associated with these sites, 6[p] = J2icp 
see figure 6 for N[p] = 2. Now, S p = n p 6, where 6 = 7r/3 is the dihedral angle of an equilateral 
triangle, and n p being an integer (including zero). I.e. 6 P obtains values in the discrete set 
{... — |7r, — |7r, 0, |7r, |7r, ...} only. Clearly, there are many smooth geometries that contain 
closed paths with associated deficit angles that are not in this set. 

Related, there is already no smooth approximation of the 2D sphere by a triangulation 
with equilateral triangles. The total number of links Ni in the lattice equals Ni = 3iV — 6, 
with iVo being the total number of vertices. On the other hand N\ = \ Y^i n ii where n 8 
denotes the coordination number of vertex i. This gives 

= 67V„ - 12, (13) 

i 

and implies restrictions for the deficit angles. For each vertex the deficit angle (i.e. the 
curvature) is simply given as Si = (6 — rii)6. Now the positive curvature of the 2-sphere at 
every point cannot be realized with equilateral triangles since n 8 < 6 is not possible for all 
vertices due to equation (13). Even if we try n 8 < 6 for each vertex, we see that (almost) 
all vertices have to have coordination number 6 and a maximum of 12 vertices might show 
positive curvature. To distribute vertices of positive curvature uniformly on the 2-sphere 
one has to face as many vertices of negative curvature. As an alternative we could collect 
all negative curvature in few (singular) vertices. 

It seems that one needs ensemble averages over equilateral triangulations to reproduce 
classical properties of geometries as simple as the 2D sphere. It is then difficult to imagine 
that fluctuation can be suppressed, for instance how to get a smooth sphere and Onsager 
exponents in the limit of no gravitational fluctuation? What are small fluctuations in the 
SDT framework? We think that even in 2D the SDT concept needs to be elaborated in much 
more details. 

10 



4 Summary and Conclusions 



Two dimensions may be a bad testing ground for quantum gravity concepts. Beautiful ex- 
act solutions exist. However, there may be reasons for their existence which obstruct their 
physical relevance. In essence, two dimensional gravity seems to be solvable, because (a) the 
classically dominant Einstein term does not fluctuate, and (b) there are (hidden) relations to 
solvable models in flat space. Therefore, we consider it unlikely that it will provide valuable 
guidance towards realistic quantum gravity problems. We have illustrated this point, by 
demonstrating consistency with no gravitational dressing, once the proper fixed geometry 
model is used as starting point to include the effects of geometry fluctuations. It seems that 
SDT relates to fixed geometry DT spin systems in a similar way as the Regge approach re- 
lates to the Onsager limit. Of course, numerical results cannot prove that the fixed geometry 
DT exponents are really identical with the KPZ exponents. But certainly our results show 
that fixed geometry DT introduces non-trivial interactions which lead out of the Onsager 
universality class. 

Acknowledgements: On of the authors (BB) would like to thank Herbert Hamber for 
drawing our attention to Ref. [7]. 
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Tables 



n 




3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 




26 


Gl 




1 


69 


49 


57 


43 


13 


5 


1 


6 


1 


1 


1 


1 


1 










1 


G2 







62 


81 


34 


26 


17 


14 


5 


5 


2 


2 


1 











1 








Table 1: Distribution of coordination numbers n for geometries Gl and G2 in case of N = 
250 sites. 



Exponent 


KPZ 


Onsager 


7/(2i/) 


0.666 


0.875 


(l-/3)/M 


0.1666 


0.4375 


1/(2*) 


0.333 


0.5 


a/{2u) 


-0.333 






Table 2: KPZ versus Onsager critical exponents. 



N 


Sweeps 


Xc 








Df 


/3 C 3 


250 


4 • 10 6 


28.60 (12) 


1.1819 (13) 


1.6290 (50) 


1.2096 (14) 


1.515 (14) 


1.0875 (32) 


500 


4 • 10 6 


46.53 (25) 


1.1720 (09) 


2.0253 (86) 


1.1968 (12) 


1.916 (23) 


1.1027 (31) 


1000 


4 • 10 6 


73.03 (51) 


1.1587 (11) 


2.400 (16) 


1.1806 (15) 


2.450 (33) 


1.1098 (19) 


1500 


4 • 10 6 


95.75 (79) 


1.1560 (11) 


2.620 (21) 


1.1773 (16) 


2.741 (49) 


1.1031 (53) 


2000 


6.7 • 10 6 


115.3 (1.0) 


1.1546 (12) 


2.832 (25) 


1.1739 (13) 


2.918 (48) 


1.1092 (36) 


3000 


10 7 


147.3 (1.3) 


1.1481 (11) 


3.039 (16) 


1.1665 (16) 


3.324 (45) 


1.1028 (26) 


4500 


7- 10 6 


194.6 (2.1) 


1.1463 (11) 


3.349 (36) 


1.1627 (19) 


3.827 (90) 


1.1045 (35) 



Table 3: Statistics and numerical results for geometry Gl. 
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Nn 


k_J VV V> L/O 


V 
AC 




^ c 


8 2 


^ c 


8 3 


250 


4 • 10 6 


28.65 (12) 


1.2138 (22) 


1.2907 (44) 


1.2094 (31) 


1.265 (15) 


1.0479 (41) 


500 


4 • 10 6 


44.24 (33) 


1.1934 (25) 


1.4491 (89) 


1.1899 (41) 


1.533 (38) 


1.004 (15) 


1000 


4 • 10 6 


67.98 (85) 


1.1769 (26) 


1.601 (13) 


1.1753 (46) 


1.736 (51) 


1.0301 (75) 


1500 


4 • 10 6 


85.4 (1.7) 


1.1722 (53) 


1.647 (19) 


1.175 (12) 


1.923 (52) 


1.0051 (71) 


2000 


6.7 • 10 6 


102.5 (1.6) 


1.1719 (41) 


1.690 (17) 


1.1770 (89) 


1.991 (35) 


1.032 (14) 


3000 


10 7 


133.2 (1.8) 


1.1670 (74) 


1.758 (15) 


1.1493 (50) 


2.104 (28) 


1.022 (12) 



Table 4: Statistics and numerical results for geometry G2. 



Observable 


Gl KPZ 


Gl Onsager 


Gl f3 c 


G2 KPZ 


G2 Onsager 


G2 f3 c 


Xc(iVo) 


0.32 


10 -25 


1.1226 (17) 


0.66 


2 • 10" 4 


1.1227 (43) 


D?(N ) 


0.52 


10" 4 


1.1319 (26) 


0.04 


7 • 10" 4 


1.1135 (73) 




0.20 


10" 3 


1.0961 (53) 


0.43 


0.20 


1.103 (63) 



Table 5: Goodness- of -fit Q , assuming KPZ versus Onsager critical exponents. In addition 
f3 c estimates are given which assume the KPZ exponent v = 1.5 



N 


Gl Cv, c 


Gl Cv(Pc) 


Gl f3 4 c 


G2 Gv,c 


G2 C v (f3 c ) 


G2 f3 4 c 


250 


0.29239 (61) 


0.22742 (65) 


1.2670 (20) 


0.25617 (59) 


0.22038 (66) 


1.2779 (49) 


500 


0.30279 (82) 


0.23449 (49) 


1.2413 (22) 


0.25427 (61) 


0.22561 (79) 


1.2569 (65) 


1000 


0.30698 (85) 


0.24291 (99) 


1.2250 (26) 


0.2521 (11) 


0.22802 (81) 


1.2618 (98) 


1500 


0.30585 (84) 


0.2458 (11) 


1.2148 (32) 


0.25337 (86) 


0.2263 (11) 


1.2563 (87) 


2000 


0.3084 (11) 


0.24650 (88) 


1.2108 (41) 


0.25212 (56) 


0.22737 (61) 


1.2325 (99) 


3000 


0.3083 (11) 


0.2499 (10) 


1.2043 (30) 


0.25219 (68) 


0.22875 (48) 


1.2239 (30) 


4500 


0.3089 (11) 


0.25017 (61) 


1.1937 (51) 








KPZ 


0.31 


0.09 




0.67 


0.014 




Onsager 


0.19 


2 • lO" 3 




0.47 


10" 3 





Table 6: Numerical results for the specific heat (geometries Gl and G2) The last two rows 
give the goodness- of -fits Q for the KPZ and the Onsager scenarios. 



Exponent 


Gl 


G2 


7/(2i/) 


6: 0.6489 (43) [0.53] 


6: 0.6185 (42) [0.48] 


(l-/3)/M 


6: 0.2211 (36) [0.45] 


5: 0.1082 (54) [0.07] 


1/(2!/) 


7: 0.3185 (52) [0.12] 


5: 0.176 (15) [0.66] 



Table 7: Two parameter fits for the exponents. The number of data used is indicated first, 
estimates and their error bars follow, finally the number [Q] gives the goodness- of -fit. 
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Figure 1: The left picture shows the flip of a link kl — > ij , preserving the 
geometry. On the right a flip is depicted which is not performed since the 
geometry changes, the area would increase. 
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1.05 1.1 1.15 1.2 1.25 1.3 



Figure 3: Geometry Gl: Susceptibility x(fi)-> derivative of the (absolute) 
magnetization D m (f3) and specific heat CV(/3). From up to down the curves 
correspond (in each case) to N = 4500, 3000, 2000, 1500, 1000, 500 and 
250. 
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Figure 4: Geometry Gl: Derivative of the Binder cumulant D B (f3). From 
up to down the data correspond to N = 4500, 3000, 2000, 1500, 1000, 500 
and 250. 
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Figure 5: KPZ versus Onsager fits for geometries Gl and Gl. Depicted is 
[Xc(-^o) ~~ Xi(-^o)]) where Xc are the measured maxima of the susceptibility, 
and xi i s the geometry Gl KPZ fit for this quantity. 
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Figure 6: A path encircles two vertices. Any vector transported around this 
path will be rotated by an angle 6 p . The trip around vertex v\ leads to a 
rotation 6 Vl , followed by an additional rotation caused by the deficit angle 8 V2 
of vertex v 2 . The total deficit angle 6 P is thus the sum 6 Vl -\-6 V2 . Stretching the 
path to encircle more vertices generalizes the result for an arbitrary number 
of vertices. 
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